Linkage disequilibrium levels in Bos indicus and Bos taurus cattle using medium and high density SNP chip data and different minor allele frequency distributions

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Strategies for single nucleotide polymorphism (SNP) genotyping to enhance genotype imputation in Gyr (Bos indicus) dairy cattle: Comparison of commercially available SNP chips

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Saliva-derived DNA performs well in large-scale, high-density single-nucleotide polymorphism microarray studies

As of June 2009, 361 genome-wide association studies (GWAS) had been referenced by the HuGE database. GWAS require DNA from many thousands of individuals, relying on suitable DNA collections. We recently performed a multiple sclerosis (MS) GWAS where a substantial component of the cases (24%) had DNA derived from saliva. Genotyping was done on the Illumina genotyping platform using the Infinium Hap370CNV DUO microarray. Additionally, we genotyped 10 individuals in duplicate using both saliva- and blood-derived DNA. The performance of blood- versus saliva-derived DNA was compared using genotyping call rate, which reflects both the quantity and quality of genotyping per sample and the “GCScore,” an Illumina genotyping quality score, which is a measure of DNA quality. We also compared genotype calls and GCScores for the 10 sample pairs. Call rates were assessed for each sample individually. For the GWAS samples, we compared data according to source of DNA and center of origin. We observed high concordance in genotyping quality and quantity between the paired samples and minimal loss of quality and quantity of DNA in the saliva samples in the large GWAS sample, with the blood samples showing greater variation between centers of origin. This large data set highlights the usefulness of saliva DNA for genotyping, especially in high-density single-nucleotide polymorphism microarray studies such as GWAS.

Promise and pitfalls of the Immunochip

Genomewide association studies (GWAS) have proven a powerful hypothesis-free method to identify common disease-associated variants. Even quite large GWAS, however, have only at best identified moderate proportions of the genetic variants contributing to disease heritability. To provide cost-effective genotyping of common and rare variants to map the remaining heritability and to fine-map established loci, the Immunochip Consortium has developed a 200,000 SNP chip that has been produced in very large numbers for a fraction of the cost of GWAS chips. This chip provides a powerful tool for immunogenetics gene mapping.

Genomic prediction for tuberculosis resistance in dairy cattle

Background: The increasing prevalence of bovine tuberculosis (bTB) in the UK and the limitations of the currently available diagnostic and control methods require the development of complementary approaches to assist in the sustainable control of the disease. One potential approach is the identification of animals that are genetically more resistant to bTB, to enable breeding of animals with enhanced resistance. This paper focuses on prediction of resistance to bTB. We explore estimation of direct genomic estimated breeding values (DGVs) for bTB resistance in UK dairy cattle, using dense SNP chip data, and test these genomic predictions for situations when disease phenotypes are not available on selection candidates. Methodology/Principal Findings: We estimated DGVs using genomic best linear unbiased prediction methodology, and assessed their predictive accuracies with a cross validation procedure and receiver operator characteristic (ROC) curves. Furthermore, these results were compared with theoretical expectations for prediction accuracy and area-under-the-ROC- curve (AUC). The dataset comprised 1151 Holstein-Friesian cows (bTB cases or controls). All individuals (592 cases and 559 controls) were genotyped for 727,252 loci (Illumina Bead Chip). The estimated observed heritability of bTB resistance was 0.23±0.06 (0.34 on the liability scale) and five-fold cross validation, replicated six times, provided a prediction accuracy of 0.33 (95% C.I.: 0.26, 0.40). ROC curves, and the resulting AUC, gave a probability of 0.58, averaged across six replicates, of correctly classifying cows as diseased or as healthy based on SNP chip genotype alone using these data. Conclusions/Significance: These results provide a first step in the investigation of the potential feasibility of genomic selection for bTB resistance using SNP data. Specifically, they demonstrate that genomic selection is possible, even in populations with no pedigree data and on animals lacking bTB phenotypes. However, a larger training population will be required to improve prediction accuracies. © 2014 Tsairidou et al.

Estudo de associação genômica e análise de enriquecimento genético da criotolerância espermática em bovinos da raça Nelore

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

The Buffalo Genome and the Application of Genomics in Animal Management and Improvement

The publication of the human genome sequence in 2001 was a major step forward in knowledge necessary to understand the variations between individuals. For farmed species, genomic sequence information will facilitate the selection of animals optimised to live, and be productive, in particular environments. The availability of cattle genome sequence has allowed the breeding industry to take the first steps towards predicting phenotypes from genotypes by estimating a genomic breeding value (gEBV) for bulls using genome-wide DNA markers. The sequencing of the buffalo genome and creation of a panel of DNA markers has created the opportunity to apply molecular selection approaches for this species.The genomes of several buffalo of different breeds were sequenced and aligned with the bovine genome, which facilitated the identification of millions of sequence variants in the buffalo genomes. Based on frequencies of variants within and among buffalo breeds, and their distribution across the genome compared with the bovine genome, 90,000 putative single nucleotide polymorphisms (SNP) were selected to create an Axiom (R) Buffalo Genotyping Array 90K. This SNP Chip was tested in buffalo populations from Italy and Brazil and found to have at least 75% high quality and polymorphic markers in these populations. The 90K SNP chip was then used to investigate the structure of buffalo populations, and to localise the variations having a major effect on milk production.

Assinaturas de seleção em bovinos da raça Canchim

Pós-graduação em Genética e Melhoramento Animal - FCAV

Imputation of non-genotyped individuals using genotyped progeny in Nellore, a Bos indicus cattle breed

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Low frequency of Y anomaly detected in Australian Brahman cow-herds

Indicine cattle have lower reproductive performance in comparison to taurine. A chromosomal anomaly characterized by the presence Y markers in females was reported and associated with infertility in cattle. The aim of this study was to investigate the occurrence of the anomaly in Brahman cows. Brahman cows (n = 929) were genotyped for a Y chromosome specific region using real time-PCR. Only six out of 929 cows had the anomaly (0.6%). The anomaly frequency was much lower in Brahman cows than in the crossbred population, in which it was first detected. It also seems that the anomaly doesn't affect pregnancy in the population. Due to the low frequency, association analyses couldn't be executed. Further, SNP signal of the pseudoautosomal boundary region of the Y chromosome was investigated using HD SNP chip. Pooled DNA of non-pregnant and pregnant cows were compared and no difference in SNP allele frequency was observed. Results suggest that the anomaly had a very low frequency in this Australian Brahman population and had no effect on reproduction. Further studies comparing pregnant cows and cows that failed to conceive should be executed after better assembly and annotation of the Y chromosome in cattle.

Genomic analysis for managing small and endangered populations: a case study in Tyrol Grey cattle

Analysis of genomic data is increasingly becoming part of the livestock industry. Therefore, the routine collection of genomic information would be an invaluable resource for effective management of breeding programs in small, endangered populations. The objective of the paper was to demonstrate how genomic data could be used to analyse (1) linkage disequlibrium (LD), LD decay and the effective population size (NeLD); (2) Inbreeding level and effective population size (NeROH) based on runs of homozygosity (ROH); (3) Prediction of genomic breeding values (GEBV) using small within-breed and genomic information from other breeds. The Tyrol Grey population was used as an example, with the goal to highlight the potential of genomic analyses for small breeds. In addition to our own results we discuss additional use of genomics to assess relatedness, admixture proportions, and inheritance of harmful variants. The example data set consisted of 218 Tyrol Grey bull genotypes, which were all available AI bulls in the population. After standard quality control restrictions 34,581 SNPs remained for the analysis. A separate quality control was applied to determine ROH levels based on Illumina GenCall and Illumina GenTrain scores, resulting into 211 bulls and 33,604 SNPs. LD was computed as the squared correlation coefficient between SNPs within a 10 mega base pair (Mb) region. ROHs were derived based on regions covering at least 4, 8, and 16 Mb, suggesting that animals had common ancestors approximately 12, 6, and 3 generations ago, respectively. The corresponding mean inbreeding coefficients (F ROH) were 4.0% for 4 Mb, 2.9% for 8 Mb and 1.6% for 16 Mb runs. With an average generation interval of 5.66 years, estimated NeROH was 125 (NeROH>16 Mb), 186 (NeROH>8 Mb) and 370 (NeROH>4 Mb) indicating strict avoidance of close inbreeding in the population. The LD was used as an alternative method to infer the population history and the Ne. The results show a continuous decrease in NeLD, to 780, 120, and 80 for 100, 10, and 5 generations ago, respectively. Genomic selection was developed for and is working well in large breeds. The same methodology was applied in Tyrol Grey cattle, using different reference populations. Contrary to the expectations, the accuracy of GEBVs with very small within breed reference populations were very high, between 0.13-0.91 and 0.12-0.63, when estimated breeding values and deregressed breeding values were used as pseudo-phenotypes, respectively. Subsequent analyses confirmed the high accuracies being a consequence of low reliabilities of pseudo-phenotypes in the validation set, thus being heavily influenced by parent averages. Multi-breed and across breed reference sets gave inconsistent and lower accuracies. Genomic information may have a crucial role in management of small breeds, even if its primary usage differs from that of large breeds. It allows to assess relatedness between individuals, trends in inbreeding and to take decisions accordingly. These decisions would be based on the real genome architecture, rather than conventional pedigree information, which can be missing or incomplete. We strongly suggest the routine genotyping of all individuals that belong to a small breed in order to facilitate the effective management of endangered livestock populations.

Prospecting major genes in dairy buffaloes

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A multiparental cross population for mapping QTL for relevant agronomic traits in durum wheat (Triticum durum Desf.)

Multiparental cross designs for mapping quantitative trait loci (QTL) in crops are efficient alternatives to conventional biparental experimental populations because they exploit a broader genetic basis and higher mapping resolution. We describe the development and deployment of a multiparental recombinant inbred line (RIL) population in durum wheat (Triticum durum Desf.) obtained by crossing four elite cultivars characterized by different traits of agronomic value. A linkage map spanning 2,663 cM and including 7,594 single nucleotide polymorphisms (SNPs) was produced by genotyping 338 RILs with a wheat-dedicated 90k SNP chip. A cluster file was developed for correct allele calling in the framework of the tetraploid durum wheat genome. Based on phenotypic data collected over four field experiments, a multi-trait quantitative trait loci (QTL) analysis was carried out for 18 traits of agronomic relevance (including yield, yield-components, morpho-physiological and seed quality traits). Across environments, a total of 63 QTL were identified and characterized in terms of the four founder haplotypes. We mapped two QTL for grain yield across environments and 23 QTL for grain yield components. A novel major QTL for number of grain per spikelet/ear was mapped on chr 2A and shown to control up to 39% of phenotypic variance in this cross. Functionally different QTL alleles, in terms of direction and size of genetic effect, were distributed among the four parents. Based on the occurrence of QTL-clusters, we characterized the breeding values (in terms of effects on yield) of most of QTL for heading and maturity as well as yield component and quality QTL. This multiparental RIL population provides the wheat community with a highly informative QTL mapping resource enabling the dissection of the genetic architecture of multiple agronomic relevant traits in durum wheat.

A deletion in the N-myc downstream regulated gene 1 (NDRG1) gene in Greyhounds with polyneuropathy

The polyneuropathy of juvenile Greyhound show dogs shows clinical similarities to the genetically heterogeneous Charcot-Marie-Tooth (CMT) disease in humans. The pedigrees containing affected dogs suggest monogenic autosomal recessive inheritance and all affected dogs trace back to a single male. Here, we studied the neuropathology of this disease and identified a candidate causative mutation. Peripheral nerve biopsies from affected dogs were examined using semi-thin histology, nerve fibre teasing and electron microscopy. A severe chronic progressive mixed polyneuropathy was observed. Seven affected and 17 related control dogs were genotyped on the 50k canine SNP chip. This allowed us to localize the causative mutation to a 19.5 Mb interval on chromosome 13 by homozygosity mapping. The NDRG1 gene is located within this interval and NDRG1 mutations have been shown to cause hereditary motor and sensory neuropathy-Lom in humans (CMT4D). Therefore, we considered NDRG1 a positional and functional candidate gene and performed mutation analysis in affected and control Greyhounds. A 10 bp deletion in canine NDRG1 exon 15 (c.1080_1089delTCGCCTGGAC) was perfectly associated with the polyneuropathy phenotype of Greyhound show dogs. The deletion causes a frame shift (p.Arg361SerfsX60) which alters several amino acids before a stop codon is encountered. A reduced level of NDRG1 transcript could be detected by RT-PCR. Western blot analysis demonstrated an absence of NDRG1 protein in peripheral nerve biopsy of an affected Greyhound. We thus have identified a candidate causative mutation for polyneuropathy in Greyhounds and identified the first genetically characterized canine CMT model which offers an opportunity to gain further insights into the pathobiology and therapy of human NDRG1 associated CMT disease. Selection against this mutation can now be used to eliminate polyneuropathy from Greyhound show dogs.

A Missense Mutation in the SERPINH1 Gene in Dachshunds with Osteogenesis Imperfecta

Osteogenesis imperfecta (OI) is a hereditary disease occurring in humans and dogs. It is characterized by extremely fragile bones and teeth. Most human and some canine OI cases are caused by mutations in the COL1A1 and COL1A2 genes encoding the subunits of collagen I. Recently, mutations in the CRTAP and LEPRE1 genes were found to cause some rare forms of human OI. Many OI cases exist where the causative mutation has not yet been found. We investigated Dachshunds with an autosomal recessive form of OI. Genotyping only five affected dogs on the 50 k canine SNP chip allowed us to localize the causative mutation to a 5.82 Mb interval on chromosome 21 by homozygosity mapping. Haplotype analysis of five additional carriers narrowed the interval further down to 4.74 Mb. The SERPINH1 gene is located within this interval and encodes an essential chaperone involved in the correct folding of the collagen triple helix. Therefore, we considered SERPINH1 a positional and functional candidate gene and performed mutation analysis in affected and control Dachshunds. A missense mutation (c.977C>T, p.L326P) located in an evolutionary conserved domain was perfectly associated with the OI phenotype. We thus have identified a candidate causative mutation for OI in Dachshunds and identified a fifth OI gene.